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We present the theory of two-dimensional, clean quantum antiferromagnets with a 
small, positive, zero temperature (T) stiffness p s , but with the ratio ksT/p s arbitrary. 
Universal scaling forms for the uniform susceptibility (xu), correlation length(£), and 
NMR relaxation rate (1/71) are proposed and computed in a 1/N expansion and 
by Monte-Carlo simulations. For large ksT/p s , Xu(T)/T and T£(T) asymptote to 
universal values, while 1/71 (T) is nearly T- independent. We find good quantitative 
agreement with experiments and some numerical studies on La^-iSr^CuO^. 
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The last few years have seen extensive theoretical and experimental studies of two di- 
mensional quantum Heisenberg antiferromagnets, with particular attention to the antiferro- 
magnetism in the cuprate compounds J5,0- On the theoretical side, most notable has been 
the work of Chakravarty et.al. who focused mainly on the low temperature (T) proper- 
ties of systems with well established long-range Neel order at T = 0; their most detailed 
results were in a regime in which the fully renormalized, T = 0, spin- stiffness p s , was not 
too small, while the temperature satisfied ksT <C p s . Under these conditions, the antiferro- 
magnet could be treated as a classical system, with all effects of quantum fluctuations being 
absorbed into renormalization of the couplings. At low T, there has been good agreement 
between their results and experiments on La^CuO^ 0. However, the experimental results 
at higher T remain poorly understood - there are clear deviations from the classical behavior 
and it is expected that quantum fluctuations will play a more fundamental role. Besides, in 
the lightly-doped cuprates, p s is likely to be quite small, thus decreasing the T range over 
which the renormalized-classical behavior will hold. Finally, there are experimental realiza- 
tions of frustrated two-dimensional Heisenberg antiferromagnets which, in all likelihood, 
have a very small value of p s . 

Our understanding of the experiments would clearly be improved by precise theoretical 
predictions in low temperature regimes other than ksT <C p s . To this end, we discuss here 
some universal properties of clean two-dimensional quantum Heisenberg antiferromagnets, 
with nearest-neighbor exchange J, in which the stiffness p s is 'small', but non-zero. We 
will study the physics when < p s < J, /c#T <C J, but the ratio /c#T '/ ' p s is allowed to be 
arbitrary. The system is then controlled by renormalization- group flows near the the T = 
quantum fixed-point separating the Neel ordered and quantum-disordered phases. Our main 
new result will be that, in this regime, the absolute values of the entire long-wavelength, 
low-frequency, uniform and staggered spin susceptibilities are completely universal functions 
of just three thermodynamic parameters: p s , c, and the ordered staggered moment Nq. The 
universal functions depend only on the symmetry of the order parameter, and sensitivity to 
all lattice-scale physics arises only through the values of p s , c, and Nq. For small ksT '/ p s , 
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the T-dependence of our results is similar to those already obtained in Ref. [[|. For large 
/cbT/ p s , most of our results are new. We will show that they are consistent with the available 
experimental 0,§],|7| and some of the numerical B|^,|ICi|j data on the uniform susceptibility, 
correlation length and NMR relaxation rate for undoped and weakly-doped Lai-bSr&CuO^. 
We thus argue that the use of a 'small' p s point-of-view is not unreasonable even for the pure 
square lattice, spin- 1/2, Heisenberg antiferromagnet; while ordered at T = 0, this system is 
evidently close to the point where long-range-order vanishes. 

Our results follow from some very general properties of the T = quantum fixed-point 
separating the magnetically-ordered and quantum-disordered phases. These properties are 
expected to be valid in both undoped and doped antiferromagnet s, though not in the pres- 
ence of randomness [fLl] , [12]1 . They are (i) the fixed point is described by a continuum 2+1 
dimensional field theory which is 'Lorentz'-invariant, and the spin-wave velocity, c, remains 
non-singular through the phase transition (ii) at T = 0, on the magnetically ordered side, 
there is a Josephson correlation length £j which diverges at the quantum fixed-point; near 
this fixed point p s equals TicQ/^j where G is a universal number []I3"|,H|; and (m) Turning 



on a small T places the critical field theory in a 'slab' geometry which is infinite in the two 
spatial directions, but of finite length L T = hc/(kBT), in the imaginary time (r) direction - 
its consequences therefore follow from finite-size scaling. 

Uniform susceptibility, \u- We first consider the response of the antiferromagnet to a 
static, spatially uniform, external magnetic field (the extension to a field at finite wavevector 
k or frequency uo will be omited here for brevity). Such a field causes a uniform precession 
of all the spins, which can be removed by transforming to a rotating reference frame at the 
price of a twist in the boundary conditions along the r direction The response of the 



system to this twisted boundary condition defines a stiffness, p T , which equals Xu- However, 
the fixed point is Lorentz invariant, and hence \ u has the same scaling properties as p s . 



Application of finite-size scaling [[RJ then yields the following T dependence for \u 



where gps/fr is the gyromagnetic ratio, N is the number of components of the order pa- 
rameter, and £Iq(x) is a universal function. Note x oc £j/L T , the length-ratio expected in 
finite-size scaling functions. We have computed Qq(x) in a 1/N expansion for the O(N) 
non-linear sigma model in 2+1 dimensions [|I7]]. The 0(3) model describes the low-energy 
dynamics of two-dimensional Heisenberg antiferromagnets on a square lattice. The antifer- 
romagnet also carries Berry phases, not present in the cr-model, but these have been argued 
to be irrelevant at the quantum fixed-point At N = oo, the scaling function Qq(x) can 
easily be calculated: 



^ 7rx 2ire l ' x 



F(x) = 2 arcsinh((l/2)e" 1/:r ) (2) 
Of particular interest is the behavior of Xu for large x. The function Qq =oq (x) is analytic 



at x = oo, and the general principles of finite-size scaling [HJ suggest that this remains true 
at finite N. Thus we expect that Qq(x — > oo) = + Qi/x + ■ ■ ■ with Q^, Qi universal 
numbers. Combined with ([I]), this implies that a plot of Xu{T) vs. T will be straight line 
at large T/p s with universal slope and intercept, whose values are related to f2oo and f^i 
respectively. At iV = oo we obtain from (g), = (y/E/n) In [(^5 + l)/2] w 0.3425 and 
Qi = 4^00/5. We have computed the first 1/N correction to and indeed found that it 
is a universal, regularization-independent number. We obtained 

tt Q (x = 00) = fioo = 0.3425(1 - 0.63/A^ + ...). (3) 

We have also performed Monte- Carlo simulations of a classical D = 3 Heisenberg ferro mag- 
net on a cubic lattice, whose phase transition is expected to be in the same universality class 
as the 0(3) sigma model. We used a lattice of size Lx LxL T (L < 30, L T < 10) at its known 
critical coupling |L8] and computed p T . It then follows that Qoo = limi r ^ 00 lim^oo L T p T 
where the order of limits is crucial. The result was Qoo = 0.25 ± 0.04, in good agreement 
with the 1/N result at N = 3. Finally, there is an analogy between f2oo and another univer- 



sal number discussed recently - the universal conductivity, oq, at the superfluid- insulator 



transition 



We turn next to small x. The N = oo result (|2|) gives the leading term Qq(x —>■ 0) = 
1/ (nx) which implies Xu(T = 0) = 2g 2 /j 2 B p s / [Ti 2 c 2 N) = (2/N)x±, where x± is the transverse 
susceptibility. This is in fact equal to the exact result expected from rotational averaging 
of an ordered quantum O(N) sigma model |§ - we have indeed found no corrections in 
the 1/N expansion at T = 0. Furthermore, it was shown in Ref. || that there are no In- 
dependent corrections to the isotropic Xu in a classical lattice rotator model. In contrast, for 
our quantum O(N) sigma model, the N = oo result contains a term linear in T at small T. 
In the 1/N expansion of this quantum model, the classical contributions appear as \ogx/N 
terms; however, as expected, the coefficient of \ogx/N in Xu was found to vanish. These 
results imply 

n Q (x^0) = l/(rrx) + T N + ... (4) 

where = 1/tt + 0(1/ N) is a universal number; some x dependence in at order 1/iV 2 
is not ruled out. 

Correlation Length, £: The scaling dimension of £ is -1, and the finite-size scaling result 
for £(T) is therefore: 

CHT) = (k B T/(hc))X Q (x), (5) 

where Xq(x) is a universal function. As for Xu, there are no non-universal factors on the 
right hand side. The numerical results for Xq depend on the precise definition chosen for 
£: we use £ 2 = — (l/S(k))(dS/dk 2 )\k=o where S(k) is the equal-time staggerd spin structure 
factor (this definition is slightly different from that used in Ref. 0). At N = oo, we found 

= F(x)y/2/[l + F(x)/sinh(F(x))] (6) 

For large x, the properties of £ -1 are similar to those of Xu- The function Xq(x) is 
expected to be analytic at x = oo with Xq(x — ► oo) = X^ + X\/x + . . .; a plot of £ _1 (T) 



vs. T will be straight line at large T with a universal slope and intercept, whose values are 
related to X^ and X x respectively. At N = oo we find X^ = 0.998, X x = —0.990. 

For small x, the N = oo result gives Xq(x — > 0) = e~ l l x . However, unlike Xu, 
the renormalized-classical spin fluctuations make a strong contribution to £, thus requiring 
careful consideration of the logx/N terms in the 1/N expansion. We identified terms to 
order (logx/N) 2 , exponentiated them and found 

X Q ( X -> 0) = ^ x -l/(AT-2) e -iV/((iV-2)s) (7) 

where Yn = is a universal constant, and liniTv^oo Yn = 1. The analysis of Chakravarty et. 
al. 0, valid for k^T <C p s but p s /J arbitrary, obtained an identical functional form for 
£(T), but with Y N replaced by a non- universal, regularization dependent prefactor whose 
determination for a square-lattice 2D S = 1/2 Heisenberg antiferromagnet required a lengthy 
calculation. Our new result here is that, in the lowest order in p s /J, this prefactor is also 
universal, and can therefore be evaluated in any convenient regularization scheme. 

NMR relaxation rate, 1/Ti: The relaxation of nuclear spins coupled to the antiferro- 
magnetic order parameter (e.g. Cu nuclear spins in La2CuO±) is given by \/T\(T) = 
linv^o 2A^ksT / (Too) J d 2 k/(4:7i 2 )x"(k,uj), where x"(k,u) is the imaginary part of the dy- 
namic staggered susceptibility of the underlying quantum antiferromagnet, and A n is the 
bare hyperfine coupling. This determines the scaling dimension of 1/Ti at the quantum 
fixed-point to be rj, the critical exponent associated with spin correlations at criticality: 
i] = 8/(3tt 2 N) — 512/(277r 4 iV 2 ) + ... in a 1/N expansion |22|] and the best current value at 



iV = 3 is f] ~ 0.028 |18| |. The finite-size scaling form for 1/Ti can be shown to be 

l/UT) = (2AlN 2 J p s )xm Q (x) (8) 



where Rq(x) is a completely universal function. Ref. [|2T| used the renormalized hyperfine 
coupling A n which is A n = A n N . Note the complete absence of non- universal normalization 
factors in (|8|). 

We now consider the limiting behavior of Rq(x) for large and small x. As before, at 
large x, Rq(x — > oo) = R^, a positive constant; the small value of r\ then implies that 1/Ti 
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is essentially T independent at high T. To leading order in 1/N, Roo can be deduced from 
the results of Ref. |TT| to be = 0.06/iV. Note the factor of 1/N - x"(k,u) is finite at 
uj — > only due to the self-energy corrections. At small x, dynamical scaling [[23[] predicts 
that Rq(x) oc £,(x). The 1/N expansion is again singular when i < 1 and we will not 
discuss it here. 

Comparison with numerical and experimental results: We have so far presented general 
scaling forms for the magnetic properties of a two-dimensional quantum antiferromagnet 
which has p s <C J. Explicit scaling functions can be calculated at iV = oo, and examination 
of 1/N corrections has been limited to those for f2oo. These corrections were however quite 
small, and we expect, in general, that 1/N expansion is robust and numerically quite accurate 
for large values of x. On the other hand, at small x, the 1/N expansion is logarithmically 
singular, and eventually changes the leading singularity in some of the scaling functions at 
x = 0; the final low-T behavior is the same as that in the renormalized-classical scaling 
theory of Ref || . The crossover between small and large x should occur for x around unity. 
Thus for x > 1 (but such that long- wavelength description is still valid), it is quite likely that 
1/N expansion will describe the experimental data better than the renormalized-classical 
theory, which, strictly speaking, is valid only for i C 1. In a square lattice, nearest-neighbor, 
S = 1/2 Heisenberg antiferromagnet, 2irp s 1.13J (Ref [p4| ) and we therefore expect that 
our large- N, large x results should work for x > 1, i.e., for k B T > 0.35J. 

The absence of any renormalized-classical corrections to Xu makes it an ideal candidate 
for testing our theory; the 1/N expansion should become accurate even at fairly small 
values of x. We start with the numerical results for Xu(T) on the square lattice S = 1/2 
antiferromagnet. There have been high-T series expansions ||, quantum Monte-Carlo || 
and finite cluster calculations [HJ. Their results all show that Xu{T) obeys a Curie- Weiss 



law at high T, reaches a maximum at k B T ~ J and then falls to a finite value at T = 
which is rather close to the rotationally averag ed 1/S result {h/gp B ) 2 Xu{T = 0) w 0.04/ Ja 2 , 
where a is the lattice spacing. For 0.35 J < T < 0.55 J, both series expansions || and Monte- 
Carlo calculations report a linear T dependence of Xu{T) (Fig .[!]). Also plotted is our 



theoretical prediction of Eqns (0,^) which, over the range of x values used in the figure, is 
well approximated by {h/gfi B ) 2 {Ja 2 Xu {T)) « 0.037a;(l + a/x), a = 0.8 + 0(1/N). This is 
remarkably close to the best fit to the data of Ref. || which gives 0.037a;(l + 0.775 fx). 

We consider next measurements of Xu(T) in weakly doped La2-sSrsCuO±. The interpre- 
tation of the experimental data even above the zero-doping T N requires caution because one 
has to subtract Van-Vleck, core and diamagnetic contributions from the measured Xu{T). 
Besides, at nonzero doping, a Pauli-like contribution from valence fermions should also be 
subtracted. Nevertheless, after subtraction was carried out, it was found [p|j25| that there is 
a large amount of universality in the susceptibility data at different doping concentrations 
and the measured susceptibility is clearly linear in T in the temperature interval 400 — 700K 
(i.e., 0.3J — 0.5J for J ~ 0.12ev). The measured slope of (h/g[iB) 2 (Ja 2 Xu(T)) vs. x is about 
0.04 which is very close to our result (0.037). 

Now the 63 Cu spin-lattice relaxation rate, 1/Ti, in La^CuO^. At very low T (x small), the 
theory of Ref. JH]] predicts that 1/Ti oc e~ 3//a; ; this is consistent with recent observations 
However, the present theory predicts that 1/Xi becomes nearly T independent for x > 1 or 
T > 0.35 J. This has in fact already been observed in series expansions || and finite cluster 
calculations JTU] for the square lattice antiferromagnet. More importantly, a flattening in 



1/T\(T) has recently been observed in the experiments on La2CuO± 0. We calculated, 
from our results above, the limiting large T value of for the same values of parameters 
as were used in the low-T fit and found 1/Xi ~ 3.3 x 10 3 sec _1 ; this is in good agreement 
with the experimental result 1/Ti ~ 2.7 x 10 3 sec _1 . Furthermore, the experimental T range 
over which \fT\ is nearly T independent increases upon doping. This is also consistent with 
our results because p s is expected to decrease with doping, thus pushing the system into 
larger x for the same T. 

The agreement between our theory and the data for £ is however not so good. Detailed 
measurements of £ in La2CuOi L have been performed at low T, where the system is in 
the renormalized-classical region ||. At the highest experimentally accessible T (= 560K 
for J = 1460K), our result = 0.023A 1 is not far from the experimental value of 
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= 0.03A 1 . We also compared our results with the numerical data for £ at higher T 
PJ26[|. For T > 0.35 J this data obeys quite well oc x(l — j/x) where 7 is close to 
one. However, the overall factor in £ in the fit is close to twice our N = 00 result. This 
discrepancy is probably due to the fact that the strong singular corrections in Xq(x) at 
small x, cause the crossover from small to large x behavior to occur dbt £l lcir ger x. Note 
however that the Monte-Carlo calculations in the quantum-critical region \TJl yield the value 
of = 1.25, which is close to our N = 00 result of X^ = 0.998. 

To conclude, we have considered in this paper the magnetic properties of two-dimensional 
quantum antiferromagnets. We argued that there exists a T range where the classical 
low-T description is no longer valid and the behavior of observables is governed by the 
renormalization-group flows near the T = quantum fixed-point. A comparison with the 
experimental data for the uniform susceptibility and e3 Cu spin-lattice relaxation rate shows 
that the intermediate behavior has been observed in the range 0.35J < T < 0.55J in 
La^CuO^. 

We thank N. Read and A. Sokol for discussions. This work has been supported by NSF 
Grant DMR88-57228. 
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FIGURES 

FIG. 1. 'Experimental' || (squares) and theoretical (line) results for the uniform susceptibility 
X u = J (ah/ 'g^B) 2 )Xu of a square lattice spin-1/2 Heisenberg antiferromagnet (a is the lattice 
spacing). There are no adjustable parameters in the theoretical result (||). Over the range of x 
plotted, the function CIq(x) is very close to its large x behavior £Iq(x) 0,^(1 + 0.8/x). We used 
this large x result with from (^) at N = 3. The theoretical and experimental slopes agree 
remarkably well. The good agreement in the intercept is somewhat surprising as its theoretical 
value (= 0.8) is known only at N = oo. 
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